Uniqueness of positive bound states to Schrodinger 
systems with critical exponents 



We prove the uniqueness for the positive solutions of the following elliptic 
systems: 



Here x £ R n , n > 3, and 1 < a,/3 < with a + (3 = 2±|. In the special case 
when n = 3 and a = 2, /3 = 3, the systems come from the stationary Schrodinger 
system with critical exponents for Bose-Einstein condensate. As a key step, we 
prove the radial symmetry of the positive solutions to the elliptic system above 
with critical exponents. 
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1 Introduction 



In this paper, we consider the uniqueness of positive solutions to the following stationary 
Schrodinger system: 

j -A(u(x)) = u(x) a v(xf m 
\ -A(v(x)) = u(xfv(x) a { > 

Here x G R n . n > 3, and 1 < a, 3 < rj± % with a + 3 = 2tt §. In the special case when 
n — 3 and a = 2, /? = 3, the systems come from the stationary Schrodinger system 
with critical exponents for Bose-Einstein condensate ([11], [15], [16], and [18]). In the 
earlier works [15], [16], and [18], people pay more attention to the elliptic system (1) 
with subcritical exponents. Very interestingly, Chen and Li have proved that the best 
constant in weighted Hardy-Littlewood-Sobolev inequality can be achieved by explicit 
radially symmetric functions (see [4] and [13]). As a consequence of their work, the 
uniqueness of positive solutions to the corresponding elliptic system (it is (1) in the 
case when a = and 3 = ^z§)) has been settled down. However, when < a, (3, the 
uniqueness of smooth positive solutions to the stationary Schrodinger system (1) is an 
open question. Generally speaking, there are very few result even for the uniqueness of 
positive solutions to the ordinary differential systems. The aim of this paper is to prove 
the radial symmetry and uniqueness of positive solutions to (1) with critical exponents 
and 1 < a < (3 < 2±|. 

As one can expect, just like in the work M.Weinstein [21] in the scalar case with 
the sub-critical exponent, that there is a closed relationship between the stationary 
Schrodinger system with critical exponent with the Hardy-Littlewood-Sobolev inequal- 
ity. As we show below, this is true. 

Since we shall use Hardy-Littlewood-Sobolev inequality to prove radial symmetry of 
our solutions, let's do an excursion about recent progress of Lieb's conjecture. Let us 
begin by recalling the well-known Hardy-Littlewood-Sobolev inequalities. Let < A < n, 
1 < s,r < oo, and ||/|| p be the L p (R n ) norm of the function /. We shall write by ||/|| p ,n 
the L p norm of the function on the domain Q. Then, classical Hardy-Littlewood-Sobolev 
inequality (HLS) states that: 



/ / ^^dxdy < C s ,,, n \\f\\r\\9\\s (2) 



for any / G L r (R n ), g G L s (R n ), and for - + - + - = 2. Hardy and Littlewood also 
introduced the double weighted inequality, which was later generalized by Stein and 
Weiss in [19] in the following form: 



f{x)9{y) dxdy 



R n JR n \x\ a ° \x — y\ x \y\P° 



< Ca ,A),s,A,n||/l|r||fi'|| S (3) 



where «o + 3q > 0, 



1-I-^<^<1-I, andl + l + A + " + /j °=2. (4) 
run r r s n 
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The best constant in the weighted inequality (3) can be obtained by maximizing the 
functional 

J(f )= [ I iMiM dxd (5) 
J bp Jr" \x\ a °\x - y\ x \y\P° X V 

under the constraints ||/|| r = \\g\\ s = 1. Then the corresponding Euler-Lagrange equa- 
tions are the system of integral equations: 



X 1 rf(x) r 1 - |^ J Rn [^pffex^y 
X 2 sg(x) s 1 = |^ J Rn \ y \4\£-y\* d y 



(6) 



where f,g>0, x G R n , and Air = X 2 s = J(f,g). 

Let u = ci/ r_1 , v = c 2 g s ~ 1 , p = q = j^i, pq ^ 1, and for a proper choice of 
constants c\ and c 2 , system (6) becomes 



u{x) = ^ J Rn j^^xdy 
v ( x ) = J^P /fl" |y| Q o| x -»| A ^ 



where u,v>0, 0<p,q<oo, < X < n, ^ < 4- < *±sa, and 4r + 4r = A+Q0+A) . 

Note that in the special case where a = and /3 = 0, system (7) reduces to the 
following system: 



„'■(,,) .... ( 8 ) 



with 

1 1 X 
+ = -. 9 

q+ 1 p+ 1 n 

It is well-known that this integral system is closely related to the system of partial 
differential equations 



f (-A)^ 2 n = v q , u> 0, in R n , 
\ (-A)t/ 2 v = vP, v> 0, in R n , 

where 7 = n — X. 

When p = q = ^3*, and u(x) = v(x), system (8) becomes the single equation: 



(10) 



n±7 

«(*)=/ T^r^dy, u>0, inR n . (11) 
jR n \x — y\ n 7 

The corresponding PDE is the well-known family of semi-linear equations 

{-Ay/ 2 u = u {n+ ^ /{n -^ , u>0, inR n (12) 
In particular, when n > 3, and 7 = 2, (12) becomes 
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-Au = u («+2)/(n-2) j u > Q> in iT. (13) 

The classification of the solutions of (13) has provided an important ingredient in the 
study of the well-known Yamabe problem and the prescribing scalar curvature problem. 
Equation (13) were studied by Gidas, Ni, and Nirenberg [8], Caffarelli, Gidas, and Spruck 
[1], Chen and Li [2], and Li [12]. They classified all the solutions. Recently, Wei and Xu 
[20] generalized this result to the solutions of more general equation (12) with 7 being 
any even numbers between and n. 

Although the systems for other real values of a, (3 between and n are of interest to 
people, we shall only concentrate in this paper to the system (1) with critical exponents 
when Ka,/?<2±|anda + /3 = 2± i. 

— ' ' — n— 2 " n— 2 

Our main results are 

Theorem 1 . Any L"^2(R n ) x L^fR") positive solution pair (u,v) to the system (1) 
with critical exponents are radial symmetric functions. 

and 

Theorem 2 . Assume that 1 < a < (3 < Then any L^(R n ) x L^(R n ) radial 

symmetric solution pair (u,v) to the system (1) with critical exponents are unique such 
that u = v . 

We point out that when u = v , the elliptic system (1) reduces to the elliptic equation 
with critical exponent (13). Then u = v is in a special family of functions: 

*** (x) = c( FTF^f )( "" 2)/2 (14) 

where t > 0, x Q G R n , with some positive constants c such that each (f> Xo ,t(x) solves (13). 
This family of functions are important in the study of (1). 

Our results are motivated from the previous work [6], where Chen, Li, and Ou con- 
sidered more general system (8) and established the symmetry and monotonicity of the 
solutions. In [3], Chen and Li also obtained a regularity result of the solutions to (8). 
To establish the symmetry of the solution to (8), Chen, Li, and Ou [5] [6] [7] introduced 
a new idea, an integral form of the method of moving planes. It is entirely different 
from the traditional method used for partial differential equations. Instead of relying on 
maximum principles, certain integral norms were estimated. The new method is a very 
powerful tool in studying qualitative properties of other integral equations and systems. 
In fact, following Chen, Li, and Ou's work, Jin and Li [9] studied the symmetry of the 
solutions to the more general system (7). 

Chen and Ma [17] discussed the Liouville type theorem for the positive solutions to 
the elliptic system (10). 

In this paper, we first prove the radial symmetry of the solutions to (1) with critical 
exponents. It is obvious that the radial symmetry of the solutions reduces (1) to a system 
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of ODEs, which has special solution pair (<f> 0jt (x) , <f> 0j t(x)) . To prove the uniqueness, we 
prove that u(0) = v(0). Then by the uniqueness of the initial value problem for ODE, 
we conclude that u = v = <f) j- This is the key observation in establishing the uniqueness 
of positive solutions for (1) with critical exponents. 

Theorems 1 and 2 will be proved in the next two sections. 



2 Proof of the Radial symmetry 

We use the moving plane method introduced by Chen-Li-Ou in [5]. We shall use the 
Hardy-Littlewood-Sobolev inequality: 

where C(n,p) is a uniform positive constant and 

Tf(x)=[ \x-y\ 2 - n f(y)dy. 

The Proof of Theorem 1. For each A G R, we denote by 

H x = {xe R n ; Xl < A}. 
For each x = (xi,:r')R n , we let 

x\ = (2A — xi, x') 

be the reflection point of x with respect to the hyperplane dH x . We let e.\ — (1, 0, ...0). 
We define 

u x (x) = u(x x ), B\ = {x G H x ; u x (x) > u(x)}, 

and 

v x (x) = v(x x ), B x = {i£ H x ; v x (x) > v(x)}. 

To do the moving plane method, we need the following formula, which is obtained 
by a change of variables. 

u(x) = / j——rdy + / , A AV f ; dy, 

Jh x \x — y\ n 1 Jh x \x x — y\ n 1 

and 

J h x \x x — y\ n ~ 2 Jhx \x — y\ n ~ 2 

Then we have 

«a(*) - u(x) = J H yA - u ^){y){^^r 2 - ^J^ V- (16) 
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Note that for x G H\, we have 



> 



\x — y\ n 2 \x\ — y\ n 2 



Then for x G B™, we have 



< u x (x) — u(x) 

< <*Ibs U \ V }^ U) dy + M Bl Ut f-^ V) dy (17) 
:=/ + // 



Let p = Using the Hardy-Littlewood-Sobolev inequality (15) we can bound the 
first term / in (17) by 

| |/| p <c(n,p)K^K-«)|^ 

I < C(n,p)|w A |^- 1 |v A |^|M A - m| p 

Here the integrations are over the set B%. 

Using again the Hardy-Littlewood-Sobolev inequality (15) we can bound the first 
term // in (17) by 

| \H\p < C(n,p)|uX _1 ( v A ~ v )\^ ^ 
I < C(n,p)\u x \p\v\\P- l \v x - v\ p 

Here the integrations are over the domain B\. Hence, we have 

K - u\ PtB « < C{n,p){\u x \p~ B \\vx\p )B u\u x - u\ P)B u + \ux\p :B v\v x \%~ B \\v\ - v\ PtB ») (20) 
Similarly, we have for following formulae for v and v x . 

V{X)= ihl d y+ I |n-2 ^ 

and 

Then we have the following estimate 

\v x - v\ p < C(n,p)(\vx\p~ B \\ux\^ B v\vx - v\ PtB « + KIp^Kfc^K - u\ p>B u) (21) 

After these preparations, we can use the moving plane method as developed in [5] to 
prove the radial symmetry of the solutions. 

At first, let's start the plane from the infinity. Indeed, for A >> 1 large enough, we 
know that the quantities 

\vx\p,B%, \U\\ Pj b%, \vx\p,B v x 
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and 

\ u \\p,B v x 

all are small, which give us that 



and 



\U\ - u\ PtB u < -\v X - v\ PtB » 



These imply that \u x — u\ p b u = and \v\ — v\ p b v = 0. These say that £?" = (/) and 
= 0. " 

Next we define 



Bl 



A = {A G R; B£ = <p for all A' > A}. 

Then it follows from the fact that u(x) — > as \x\ — > oo and u(x) > in R n that 
A < +oo. By the definition of A , we have u Xo (x) < u(x) foe x G H\ . Using the 
expression (16), we see that u\ () (x) < u(x) for x G H\ . This implies that |2Aei — B™\ — > 
as A — > A . It is now standard to know (see [5]) that ii\ = u, which then gives us v\ = v. 
Since X\ can be any directions, we conclude that u and v are radial symmetric about 
some point x a . 



3 Proof of the Uniqueness 

In some sense, the proof of Theorem 2 is just at hand by using the integral expression 
of the solution pair (u, v). 

2n 2n 

Proof of Theorem 2. Let (u,v) G L™- 2 (R n ) x L n - 2 (R n ) be a pair of solutions to 
system (1). By Theorem 1, we know that u and v are radial symmetric about the some 
point xq. We may say, x = 0. 

Recall that we have assumed that 1 < a < (3 < Since, u G L~(R n ) and 

v £ L^2(R n ). Using the same method in [5], we have that u G C 2 (R n ) and v G C 2 (R n ) 
with 

u(x) — > 0, v(x) — > 0, 

as |x| — > oo. 

Since our solution u is radially symmetric, hence we can write, in polar coordinates, 
the first equation in (1) as 

(r n -V(r))' = -r n -yr) c \>(r) /3 , 

where r = \x\. 

Integrating both sides from to r yields 

r n -V(r) = - / " s n - 1 u a v (3 (s)ds. 
Jo 
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It follows by another integration that 

u(r) = u(0) - I" ~~r ( T s n - l u a v p dsdr. (22) 

JO T n JO 

Similarly, for v(r), we have 

v ( r ) = v(0) - [ T ~~r f T s n - 1 v a u p dsdr. (23) 

JO T n JO 

As we mentioned in the introduction, we need only to show that u(0) = v(0). Oth- 
erwise, suppose 

u(0) < v(0), (24) 
then by continuity, for all small r > 0, 

u(r) < v(r). (25) 

In other word, there exists an R > 0, such that 

u(r) < v(r), Vr G (0,R). (26) 

Let R Q be the supreme value of R, such that (26) holds. Then R a < oo and u(R ) = 
v(R ), where we have used the fact that u(+oo) = v(+oo) = 0. By the definition of R a 
and a < (3, we have that 

u(r) a v(rf > v(r) a u(rf, Vr e (0,i? o ). (27) 

Then we have from (22) and (23) that 

> u(0) - v(0) = f Ro [ T s n -\u a vP - u p v a ){s)dsdT > 0. 

JO T n 1 JO 

This is impossible. 

Similarly, one can show that u(0) > v(0) is impossible. Therefore, we must have 

u(0) = v(0). 

Finally, by the standard ODE theory, we arrive at 

u(r) = v(r). 

Hence, our elliptic system (1) has been reduced to the elliptic equation with critical 
exponent (13). By now, it is standard to know that our solutions pair u and v are of 
the form (14). This completes the proof of Theorem 2. 
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